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The 7p — >■ TT^TT^p reaction has been measured from threshold to 1.4 GeV using the Crystal Ball 
and TAPS photon spectrometers together with the photon tagging facility at the Mainz Microtron. 
The experimental results include total and differential cross sections as well as specific angular 
distributions, which were used to extract partial-wave amplitudes. In particular, the energy region 
below the _Di3(1520) resonance was studied. 

PACS numbers: 25.20.Lj, 13.60.Lc, 14.20. Gk 



I. INTRODUCTION 

Although studied for a long time, the properties of 
many baryon resonances are still not well known and a 
clear understanding of resonances in QCD is still not pos- 
sible. Some states below 2.5 GeV are believed to couple 
strongly to final states with two pseudoscalar mesons. 
Therefore, the investigation of tttt and ttt; photoproduc- 
tion provides important new information about the nu- 
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cleon excitation spectrum. 

During the last two decades, an extensive study of 
double-pion photoproduction for i?-y < 1 GeV has been 
undertaken |l|-[l4|. The theoretical interpretation of the 
data in various isospin channels was carried out using dif- 
ferent phenomenological analyses [T5l - [20j . As a rule, the 
models for double-pion production are based on isobar 
models or effective field theories. Typically, the reac- 
tion amplitude is constructed as a sum of background 
and resonance contributions. The background part con- 
tains nucleon Born terms as well as meson exchange in 
the t channel. The resonance part is a coherent sum of 
s-channel resonances decaying into TTTriV via intermedi- 
ate formation of meson-nucleon and meson-meson states 
( "isobars" ). Despite significant qualitative differences be- 
tween the models, in general they provide an acceptable 
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description of the existing cross-section data. Such an ap- 
parent consistency between theoretical models does not 
indicate a high level in understanding double-pion pho- 
toproduction; rather it demonstrates a weak sensitivity 
of the existing data to the underlying dynamics. 

The reaction 7p — > n'^n'^p is a typical example where 
the data on double-pion production have not been fully 
understood theoretically, especially below the second res- 
onance region. A widely accepted property of this re- 
action is a large contribution from Di 3(1520), which is 
known to couple strongly to the tt A channel [2l|, [22| • The 
Z3i3(1520) contribution to the jp -> 7r°7r°p total cross 
section is seen as the first peak at Ej « 730 MeV (see 
Fig.[T]), the features of which are reproduced more or less 
successfully by all models. However, the dynamics un- 
derlying this reaction in the region from the Z3i3(1520) 
resonance down to threshold have not been well under- 
stood so far. In this region, the total cross section demon- 
strates an almost linear rise, hinting at s-wave dominance 
in the final state, which, however, is not confirmed by 
the theory. An attempt to describe such behavior by 
a large contribution of the Roper resonance Pii(1440), 
decaying into aN in s wave Il6l. seems to be ruled out 
by subsequent investigations |12l. [isl [TtI . [l9l . [20j . The 
£)i3(1520) contribution itself, according to the results of 
Refs. [T2I. [itI [19^ . reduces rapidly with decreasing energy 
and cannot explain the experimental data in the region 
below = 650 MeV. In Ref. 12], the authors try to 
describe the jp — tt'^tt'^p reaction by a dominant contri- 
bution from the 1533(1700) resonance. The well-known 
minimum at W = 1.6 GeV and the second maximum at 
W = 1.7 GeV, seen in the 7p — ^ ir^ir^p total cross section, 
were described in Ref. [IJ] by the interference between 
Di3(1520) and 1)33(1700). However, a simple consider- 
ation using the Clebsch-Gordan coefficients shows that 
ttA photoproduction in the 7 = 3/2 channel should lead 
to the ratio 

a{7r+n~p) _ cr(7r~ A++) + cr(7r+A°) 
CT(7r"7rV) CT(7r"A+) 

{/ 31 33 \2 / 31 3_l \ 2 -| 

/ 3 i 3 1 \ 2 

: (C'ro|iCro|i j = 13, (1) 

where for simplicity we neglect the interference between 
two possible nN pairs. Taking a value a — Afih for the 
^^33(1700) contribution to n^n^p, as predicted in [l2|, as- 
sumes that at least 70% of the jp — > tt+tt^p total cross 
section comes from this resonance alone, which seems un- 
likely. The major (about 80%) part of the 7p — )■ 7r+7r~p 
total cross section comes from the A-Kroll-Ruderman 
term, so that the addition of such a strong contribution 
from 1)33(1700) would result in a significant overestima- 
tion of the experimental data. Thus the dynamics of 
double-TT*^ photoproduction in the energy region below 
1)13(1520) are still far from being well understood. 

Major disagreements between the results of different 
models as well as between the theoretical predictions and 



the experimental data were revealed in Refs. [9| and [13|, 
in which the measurement of the beam helicity asym- 
metry 7® for 7r+7r~, 7r+7r°, and 7r°7r° was reported. As 
discussed in Ref. [23], this quantity is very sensitive to the 
model details, so that even a small variation of the model 
parameters can change the results significantly. How- 
ever, the interpretation of polarization measurements in 
terms of the spin and parity of is quite difficult, es- 
pecially for the processes with more than two particles 
in the final state. Therefore, it is desirable to find a 
method that, on the one hand, will be sensitive to the de- 
tails of the dynamical structure and, on the other hand, 
will provide a clear interpretation of the results in terms 
of spin-parity of the contributing waves. Furthermore, 
especially important is that the method should not be 
connected strictly to the isobar model, allowing one to 
perform the partial-wave analysis with a minimal model 
dependence. An approach that seems to obey the re- 
quirements discussed above was applied for analysis of 
inelastic pion -nucleon scattering ttN — > ttttN (see, for 
example, |24l. [2H ) . A similar formalism for photoproduc- 
tion of two pseudoscalars was developed in Ref. [26] . Such 
approaches require high-statistics data covering the full 
solid angle. In the present work, we remeasured double- 
7r° photoproduction off the proton with an unprecedented 
accuracy and applied the formalism of Ref. 26] to study 
its dynamics, having a main goal to learn which 
waves dominate in this reaction at the energies below 
i:'i3(1520). 

The experimental data in the present study were ob- 
tained at the Mainz tagger photon facility using an al- 
most Att detector based on the Crystal Ball and TAPS 
multiphoton spectrometers. 

This paper includes a brief description of the experi- 
mental setup, data handling, the model formalism, dis- 
cussion of the results, and conclusions. 



II. EXPERIMENTAL SETUP 

The reaction jp — > tt^tt^p was studied using the Crys- 
tal Ball (CB) [23] as the central spectrometer and TAPS 
US [2^ as a forward spectrometer. These detectors 
were installed in the energy-tagged bremsstrahl ung pho- 
ton beam of the Mainz Microtron (MAMI) [13, IMl . The 
photon energies were determined by the Glasgow tagging 
spectrometer [32| - [3^ . 

The CB detector is a sphere consisting of 672 optically 
insulated Nal(Tl) crystals, shaped as truncated triangu- 
lar pyramids, which point toward the center of the sphere. 
The crystals are arranged in two hemispheres that cover 
93% of Att sr, sitting outside a central spherical cavity 
with a radius of 25 cm, which is designed to hold the 
target and inner detectors. In this experiment, TAPS 
was arranged in a plane consisting of 384 BaF2 coun- 
ters of hexagonal cross section. It was installed 1.5 m 
downstream of the CB center covering the full azimuthal 
range for polar angles from 1° to 20°. More details on 
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the energy and angular resolution of the CB and TAPS 
are given in Refs. |35l. l36j. 

The present measurement used 855-MeV and 1508- 
MeV electron beams from the upgraded Mainz Mi- 
crotron, MAMI-C Jlj. The data with the 1508-MeV 
beam were taken in 2007, and with the 855-MeV beam 
in 2008. Bremsstrahlung photons, produced by the 1508- 
MeV electrons in a 10-/im Cu radiator and collimated by 
a 4-mm-diameter Pb collimator, were incident on a 5- 
cm-long liquid hydrogen (IH2) target located in the cen- 
ter of the CB. The energies of the incident photons were 
measured in the range 617 to 1402 MeV by detecting the 
post-bremsstrahlung electrons in the Glasgow tagger 3^ . 
With the 855-MeV electron beam, bremsstrahlung pho- 
tons were produced in a diamond radiator, collimated by 
a 3-mm-diameter Pb collimator, and incident on a 10- 
cm-long IH2 target. In this experiment, the energies of 
the incident photons were tagged from 84 to 796 MeV. 
The energy resolution of the tagged photons is mostly 
defined by the width of the tagger focal plane detectors, 
and by the electron beam energy. For a beam energy 
of 1508 MeV, a typical width of a tagger channel was 
about 4 MeV, and about 2 MeV for a beam energy of 
855 MeV. Due to the beam collimation only part of the 
bremsstrahlung photon flux reached the IH2 target. In 
order to evaluate the reaction cross sections, the proba- 
bility of bremsstrahlung photons reaching the target (the 
so-called tagging efficiency) was measured for each tagger 
channel. The typical tagging efficiency in the experiment 
with the 1508-MeV electron beam was found to vary be- 
tween 67% and 71%. With the 855-MeV electron beam, 
the tagging efficiency varied with photon energy between 
30% and 60%. 

The experimental trigger in the measurement with the 
1508-MeV electron beam required the total-energy de- 
posit in the CB to exceed ^ 320 MeV and the number of 
so-called hardware clusters in the CB to be larger than 
two. With the 855-MeV electron beam, the trigger re- 
quired the total energy in the CB to exceed ~ 100 MeV, 
and the number of hardware clusters in the CB and TAPS 
together to be larger than 1, with at least one hardware 
cluster in the CB. 

More details on the experimental conditions of the data 
taking with the 1508-MeV electron beam in 2007 are 
given in Refs. [35l.[36t. 



III. DATA ANALYSIS 

The reaction 7p — tt'^tt'^p was identified using events 
with four photons detected in the calorimeters. There 
were two independent analyses made to crosscheck the 
results. In the first analysis, the event-selection proce- 
dure was similar to the one that was used to measure 
the reaction 
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37| . The second analysis was 



based on the kinematic-fit technique and was similar to 
those published in Refs. [s^, [s^. Both analyses are in 
excellent agreement. Since the kinematic- fit technique 
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1: (Color online) Total cross sections for jp — > tt^tt^p 



are shown as a function of the incident-photon energy. The 
results obtained in this work are compared to the existing data 
from GRAAL 6], CB-ELSA 0,(3, DAPHNE >], TAPS [Q, 
and Crystal Ball/TAPS [43|. Only statistical uncertainties are 
shown for all data. 



typically yields data with better resolution, the results of 
the second analysis were used. 

The details of the kinematic-fit parametrization of 
the detector information and resolution were given in 
Ref. [3^. The four- and five-cluster events that satis- 
fied the hypothesis of the process jp — 7r°7r°p — A-fp 
at the 2% confidence level, CL, (i.e., with a probabil- 
ity of misinterpretation less than 2%) were accepted as 
the reaction candidates. The kinematic-fit output for 
which the pairing combination of the four photons to 
two TT^s had the largest CL was used to reconstruct the 
reaction kinematics. Possible background from other re- 
actions was checked by their simulation and by determi- 
nation of a probability for them to be misidentified as 
27r° events. Below the 7p — >■ n^rjp threshold, the back- 
ground contamination was found to be only from interac- 
tion of incident photons in the windows of the target cell 
and from random coincidences. This contamination was 
subtracted from the experimental spectra by using data 
samples with random coincidences and with an empty 
(no liquid hydrogen) target. The background from the 
7P — n'^'qp — > 47P events was estimated to be quite 
small, reaching only 0.5% at incident-photon energies of 
1.4 GeV. So this background was neglected in our results. 

The determination of the experimental acceptance was 
based on a Monte Carlo (MC) simulation of the jp — 
tt^tt'^P reaction with different event generators based on 
various assumptions about the reaction dynamics. For 
the most part, the MC simulation was made as the pro- 
cess 7p — > Att^ — >■ 7r°7r°p, using the mass (1210 McV) 
and width (100 MeV) of the Delta resonance at its pole 
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FIG. 2: Experimental m^{Tv^''p) invariant-mass distributions 
(crosses) compared to those obtained from the MC simula- 
tion (solid line) of the 7p — >■ tt^tt^p reaction. The MC event 
generation included the Att" and -Di3(1520)7r° intermediate 
states in double- vr" photoproduction. 



position (Ref. [38|). With these parameters, the agree- 
ment between the experimental and MC-simulation dis- 
tributions of the invariant mass m('K^p) is much better 
than that obtained when the Breit-Wigner parameters 
of A from Ref. [s^ are used. The same parameters also 
give a good description of the A peak seen in the re- 
action ■n~p — )• TT^TT^n [11]. One simulation was made 
by generating an isotropic angular distribution of the 
A — > n^p decay. Another was made similar to the ex- 
perimental angular distribution in the region of the first 
peak in the — >■ 7r°7r°p total cross section (see Fig.[T|). 
Part of the MC simulation modelled the process ^p — > 
L'i3(1520)7r" — TT^TT^p, where the D13 mass and width 
were taken as 1510 MeV and 110 MeV, respectively. This 
simulation was used only in the analysis of the data taken 
with the 1508-MeV electron beam and only for the ener- 
gies in which the 7p — >■ Di3(1520)7r° — >■ ir^n^p contribu- 
tion becomes visible. For these energies, the determina- 
tion of the experimental acceptance was done by mixing 
the A and -D13 simulations, where the weights of each 
simulation were adjusted to get the best agreement with 
the experimental raiii^p) distributions. The comparison 
of the measured and simulated m^iir^p) distributions is 
shown in Fig. [2] for four different energies. The agree- 
ment between the measured and simulated distributions 
is better at lower energies. This agreement is almost in- 
dependent of the choice of the angular distribution used 
for theA-decay simulation; it affects only the shape of 
the m(7r''7r°) spectrum. The resonance peaks seen in the 



m?{'K^p) spectra look different from the Breght-Wigner 
shape as every event is represented by two m?{-n'^p) val- 
ues, which are located symmetrically in the tt^tt^p Dalitz 
plot with respect to its symmetry line (see Ref. [39] and 
its Fig. 6 for more details). Then every resonance band 
in the Dalits plot has its reflection with respect this sym- 
metry line. In Fig.[2ja) for example, the projection of 
the A band to the m?{'K^p) axis is seen in the right part 
of the spectrum, while a bump in the left part just cor- 
responds to the reflection of the A band with respect to 
the symmetry line. 

For the data at each electron-beam energy, the 
corresponding MC events were propagated through a 
GEANT (version 3.21) simulation of the experimental 
setup, folded with resolutions of the detectors and con- 
ditions of the trigger. The resulting simulated data were 
then analyzed in the same way as the experimental data. 
The average acceptance for the data with the 855-MeV 
electron beam was found to be close to 60% for the entire 
energy range of double- 7r° photoproduction. The aver- 
age acceptance for the data with the 1508-MeV electron 
beam decreases smoothly from 55% at — 617 MeV to 
42% at = 1400 MeV'. 

The total cross sections obtained from the 7p — ?> n^n^p 
reaction are shown in Fig. [T] as a function of the incident- 
photon energy and are compared to some previous mea- 
surements. The majority of previou s — ^ tt'^tt'^p exper- 
iments were at MAMI [l|-i,[3,[l[llli3- In Fig lll we in- 
clude only the most recent results obtained with three dif- 
ferent experimental setups: DAPHNE ^8], TAPS [12,], and 



Crystal Ball/TAPS W 
performed at ELS A [11 



. The other measurements were 
, [13 and by GRAAL [^. The re- 
sults obtained in this work are in good agreement with all 
previous measurements within the given statistical and 
systematic uncertainties. It was possible with our new 
data to reduce considerably the energy binning as well as 
the statistical uncertainties. The agreement of our total- 
cross-section results from the two measurements with dif- 
ferent electron-beam energies can be seen in the overlap- 
ping range from E-y = 617 MeV to E^ = 796 MeV (shown 
later in Fig. [7]). 

The systematic uncertainties in the total and differen- 
tial cross sections were estimated to be not larger than 
6% and are dominated by the determination of the ex- 
perimental acceptance for ^p tt^tt^p and the photon- 
beam flux. The systematic uncertainty because of the 
acceptance determination was studied by comparing our 
results for the total cross sections that were obtained 
with various MC simulations based on event generators 
with different 7p — ?> n^ir^p dynamics. Also, we compared 
the total cross sections that were obtained from the in- 
tegration of the differential cross sections, which will be 
shown later in the text. The systematic uncertainty in 
the photon-beam flux was determined mostly by the vari- 
ation of the tagging efficiency during the data-taking pe- 
riod. 
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FIG. 3: Definition of tfie coordinate system used in tlie 
present formalism, k, p, q{, and qi are respectively three- 
momenta of the incident-photon, out-going proton, two pions 
in the center of mass system. Axis Z is a. normal to the decay 
plane. Axis X is along p. O and $ are respectively the polar 
and azimuthal angles of k. 



IV. THE MODEL 

The formalism used to interpret our experimental data 
is described in Ref. [1^, where the formal expressions 
are derived for the helicity amplitude as well as for the 
cross section. At a particular photon energy, the reaction 
amplitude is determined by four independent continuous 
variables for unpolarized experiments, described below, 
and two discrete variables, which are taken as the initial 
and final helicities of the nucleon. Our choice of coor- 
dinate system is shown in Fig. [31 where all particles are 
in the center of mass frame. Axis Z is chosen along the 
normal to the decay plane, which is defined by the three 
final-state particles. Axis X is chosen along the outgoing- 
nucleon momentum. Angles $ and Q defined in Fig. [3] 
specify the direction of the incident-photon momentum 
k in this coordinate system. Together with the two an- 
gles, the energies of the two pions, uji and ui2, uniquely 
determine the final-state kinematics. The angular depen- 
dence of the matrix element is then given |26{ by 

JM 

X Dij^i-f,e,~-f), (2) 

where D^_^^^ are the Wigner functions and JM denote 
respectively the total angular momentum and its projec- 
tion on axis Z. The complex partial amplitudes t'ljf , 
which depend on the energies uji and i02, contain the full 
dynamics of the process. 

After spin summation and appropriate integration over 
LUi and UJ2, one obtains the unpolarized differential cross 



section da /(d cos Qd^) or the corresponding normalized 
quantity |26| 



T¥(e,$) 



1 da 
a dQ 



E E 

L>0 M=-L 



2J+1 

An 



WLAiYLMie,^), (3) 



which is expanded over spherical harmonics with M^oo = 
1. The coefficients Wlm in Eq. (jS)) are hermitian combi- 
nations of the partial- wave amplitudes t'l^^ . The corre- 
sponding expression was obtained in [26] : 

Wlm = ^IC f dio,dL02j2 E (-1)*'+'' 

JJ'MjM'j 

where K, is an appropriate phase space factor. Formula 
determines the general structure of an angular dis- 
tribution in a manner analogous to the expansion of the 
cross section for single-meson photoproduction in terms 
of the Legendre polynomials. 

To limit the number of model parameters, only the 
lowest partial waves were used. Their choice is motivated 
by previous isobar-model analyses which demonstrated 
that only waves with J < 3/2 were important below = 
0.8 GeV [H [13, 111- 

Expansion ([3]) written for the case of J < 3/2 is 




+5(vK2oP2°(cose) 



3 

1 

'71 



VFiiPii(cose)cos$ 
VF22P2 (cose)cos2$ 
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2V5 



\/3VF3iP3(cose)cos$ 



VF33P|(cOSe) COS 3$ 



(5) 



from which it is easy to see that the cos 8 distribution 
has a general form 



VK(cos Q)^A + B cos^ e , 



with 



A=-(l- -W2U 
2 V 2 " 



B = ^W^2o . 



(6) 



(7) 



As shown in Ref. |26i] , the identity of the two pions 
together with parity conservation results in the following 
symmetry relations: 



(8) 



Then, using the known properties of the spherical har- 
monics, 
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FIG. 4: (Color online) Distribution W(e) = /VK(e,$)d$ 
shown as a function of cos Q, where is the polar angle 
of the incident photon in the coordinate frame presented in 
Fig-El Our experimental results with statistical uncertainties 
are shown by filled circles. The predictions from the model 
of Ref. [m are shown by dashed lines. The results of fitting 
our data below E.y = 0.8 GeV are shown by solid lines. The 
energy label in each panel indicates the central photon energy 
for each bin. 



FIG. 5: (Color online) Distribution W(<E>) = 

TT J VF(0, $) sin 0cO, where $ is the azimuthal angle of 

the incident photon in the coordinate frame presented in 
Fig.m Other notations are the same as in Fig.|4l 



YLMie, 27r - $) - <f ) , (10) 

one can see that the relations ([8|) lead to the following 
restrictions for the coefficients Wlm 

Wlm = , if L + M = odd , (11) 

and 

Im{ Wlm) = . (12) 

In particular, Wlq = for L = 2n + 1. Furthermore, 
parity conservation requires that the amplitudes with the 



same parity interfere only in Wlm with even L, whereas 
the waves having the opposite parity interfere only in 
Wlm with odd L. This property was effectively used in 
partial wave analyses of inelastic pion-nucleon scattering 

The rule pT|) requires that, for example, the states 
with ~ i produced via All absorption (which in 
our case is saturated by the Roper resonance and by the 
major part of the Born terms) can contribute only to 
Woo- Therefore, in the region where states with J > 
3/2 are not important, the angular distribution W{Q, $) 
should be isotropic in both Q and $. 
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FIG. 6: (Color online) Moments Wlm (normalized such that 
Woo = 1) as a function of the incident-photon energy. Our 
experimental results for the real part of Wlm are shown by 
filled circles. The fit results are shown by solid lines. 



V. DISCUSSION OF THE RESULTS 

The measured distributions of the angles Q and $ de- 
fined in Sec. IIVI are shown in Figs.H] and [5] for 50-MeV- 
wide energy bins. The predictions that were made for 
these angular distributions by the model from Ref. [l^ 
are shown in the same figures by dashed lines. This model 
describes roughly the experimental distributions in the 
energy region close to -Di3(1520). At the energies below 
= 900 MeV, the shape of the measured cos Q distri- 
butions shows good agreement with formula How- 
ever, at the energies below Z3i3(1520), where the model 
of Ref. [1^ involves the Pi 1 (1440) resonance and the nu- 
cleon Born terms come into play, the model predicts an 
angular dependence that is weaker compared to the ex- 
perimental data. As follows from Eq. (jS]), the weakening 
of the moment W20 leads to the model failure at these 
energies. Within our approximation J < 3/2, this mo- 
ment is saturated by the waves corresponding to the to- 
tal angular momentum J — 3/2. Therefore, this ob- 
servation indicates the persistence of such waves (and 
perhaps higher waves) at these energies. As already dis- 
cussed in Sec. I, the Roper resonance cannot dominate at 
these energies. At the same time, the calculations from 
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FIG. 7: (Color online) Total cross section for jp tt^tt^p as 
a function of the incident-photon energy. Our experimental 
results are shown by triangles and circles, respectively for the 
data with the 855-MeV and 15Q8-MeV electron beam. Only 
statistical uncertainties are shown. The fit results for the 
total cross section are shown by the solid line, and for the 
3/2~, 3/2^, and 1/2"*" waves by long-dashed, dash-dotted, 
and dotted lines, respectively. The Di3(152Q) contribution, 
calculated from the model of Ref. [l^ , is shown by the short- 
dashed line. 



Refs. [H, [13, m predict a rapid fall of the Di3(1520) 
contribution at lower energies. This results in the signifi- 
cant underestimation of the measured total cross section, 
which demonstrates almost linear ener gy d ependence in 
this region. According to the fit of Ref. [l^ , such behav- 
ior of the experimental data was explained by a contri- 
bution from the A-like resonance Z?33(1700), which dom- 
inated the tt'^tt'^ channel in the full energy region consid- 
ered. 

Using Eq. (jS]), the expansion coefficients Wlm were ob- 
tained from the experimental two-dimensional plots of 
cos 8 versus To illustrate the partial- wave content 
of the jp — >■ tt'^tt'^p amplitude in more detail, we show 
in Fig[B] the variation of these coefficients in the energy 
range = 400 - 1400 MeV for the waves with J < 5/2. 
As remarked above, the values of the coefficients with odd 
J are determined by the interference of the states with 
different parities. If the insignificance of the waves with 
J > 3/2 is assumed, nonzero W3M coefficients arise from 
the interference between 3/2^ and 3/2+. As one can see 
in Fig.m W31 and W33 are quite small. This observation, 
for example, may point to a predominantly background 
nature of the partial wave with = 3/2+, which thus 
has a small imaginary part, whereas — 3/2^ is mostly 
imaginary because of the closeness to the Di3(1520) pole. 
Furthermore, as will be shown later, the weakness of the 
interference between 3/2 ~ and the positive-parity states 
1/2+ and 3/2+ results in a small forward-backward asym- 
metry in the angular distributions for the final-state pi- 
ous. 




a „ [ 825 MeV 



925 MeV 



q [ 975 MeV I q I 1025, MeV ] 1075, MeV 



q \ 1125MeV I q I 1175,MeV 1 q[ 1225 MeV 



0.5. 



0.5 



q [ 1275, MeV I q I 1325, MeV | ,1375 MeV 



0.5 



-1 1-1 

COS0_ 



1 -1 



FIG. 8: (Color online) 7p — >■ ir'^Tv^'p differential cross sections 
as a function of the production angle of the outgoing tv'^ in the 
center of mass frame. Since there are two identical pions, each 
cross section represents the average of two distributions. Our 
experimental results with statistical uncertainties are shown 
by filled circles. The predictions from our model are shown by 
solid lines. The dashed lines result from the Bonn-Gatchina 
model fli],[il. 



Already at low energies, the quantities W20 and W22, 
which (within our restriction J < 3/2) are determined 
exclusively by the incoherent sum of the states 3/2~ and 
3/2+, achieve relatively large values. This observation in- 
dicates an additional strong 3/2~ contribution, interfer- 
ing with £'13(1520), and/or a large fraction of 3/2+. The 
latter can come, for example, from A decaying to ttA, 
followed by A -> TriV. We cannot also exclude a strong 
^33(1700) amplitude, as was found in Ref. [l^]- How- 
ever, as discussed in Sect. I, the experimental data on 
7r+7r~ photoproduction seem to leave no room for strong 
coupling to this resonance. 

The coefficient Wn, coming from the interference of 
the 3/2~ wave with the positive-parity waves 1/2+ and 
3/2+, demonstrates quite sharp energy dependence in the 




FIG. 9: (Color online) Same as Fig.|8] but for the outgoing 
proton. 



region Ej — 500 — 650 MeV. The moments with L = 5 
are small. 

To fit the measured values of the moments Wlm, we 
assumed the model in which the final TTTriV state is pro- 
duced exclusively via the intermediate ttA state. First, 
the resonance Z3i3(1520), whose role in this reaction is 
more or less firmly established, was put into the 3/2~ 
wave. The corresponding amplitude was parametrized 
in the Breit-Wigner form, with parameters taken from 
PDG [38] . The only other partial waves included were 
those that lead to s and p wave in the final ttA state: 
= 1/2+, 3/2", and 3/2+. From our fit, the wave 
— 5/2+, containing ttA in a p state, is negligibly 
small and was excluded from further consideration. Each 
partial-wave amplitude was parametrized in the form 



tB+tRiW) GaFa^^ 



tN ■ 



(13) 



where the two terms in the brackets stand for a smooth 
background and a rapidly varying part tji{W), which 
can contain s-channel resonances. The factors Ga and 
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FIG. 10: (Color online) jp — ^ tt^tt^p differential cross sections 
as a function tlie invariant mass squared m'^{Tv°Tr^). Our ex- 
perimental results with statistical uncertainties are shown by 
filled circles. The predictions from our model are shown by 
solid lines. The dashed lines result from the Bonn-Gatchina 
model 0,[i3- 



^a-s-ttTV are respectively the propagator and the ttN- 
decay vertex of the A isobar. The energy-independent 
background in each partial wave was parametrized as 



tB = Ksle*"*- , 



(14) 



with adjustable constants \tB\ and (j)B, whereas the 
rapidly-varying part was taken as 



tR{W) = \tRiW)\e''f"^^'^^ 



with 



\tR{W)\^ao + aiq + a2q^ , (I>r{W) 



hi + b2q^ 



(15) 



(16) 



FIG. 11: (Color onhne) Same as Fig.[TOl but for m\iT°p). 
Since there are two identical pions, each cross section repre- 
sents the average of two distributions. 
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FIG. 12: (Color online) The partial wave contributions to the 
■KIT spectrum for — 425 (a) and 775 MeV (b) (enlarged 
from Fig. I10[) . Long-dashed, dash-dotted, and dotted lines 
correspond to the 3/2~, 3/2"^, and 1/2"'' waves, respectively. 
The sum of all waves is shown by solid lines. 
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where q is the maximum pion momentum in the final 
TTTriV state, corresponding to the total energy W in the 
center of mass system 



The coefficients at and bi in Eq. ([T6l) should be deter- 
mined from the fit. Since our fit was restricted to a lim- 
ited energy range, the parametrization of |tjj(W)| by a 
simple polynomial formula (1161) was expected to be sat- 
isfactory. 

The results of the fit to the Wlm moments are shown 
by solid lines in Fig.|6l Instead of listing the results 
for the fit parameters, the integrated partial cross sec- 
tions ajp are shown in Fig. [71 As expected, the strongest 
^0^0 production into the entire energy region comes from 
— 3/2 ~ and 3/2+ (shown by the long-dashed and 
dash-dotted lines in Fig. [7]). At lower photon energy en- 
ergies the partial cross section (T3/2- falls off slower than 
predicted by the model [ll]. In addition, the wave 3/2+ 
turns out to be very important, especially at the energies 
below = 650 MeV. 

After fitting the model parameters to the measured 
total cross sections and moments Wlm, the reliability 
of our parametrization [given by Eqs. (|13p -(fT6 l) ] was 
checked by comparing the model predictions with the 
7P — > tt'^tt'^p experimental results for other observables. 
These comparisons are shown in Figs. [8HT31 Before dis- 
cussing the agreement between the experimental data 
and the calculation, it is important to note that the val- 
ues of Wlm do not determine final-state distributions 
of 7p — > TT^TT^p [because of the integration over ener- 
gies oJi and UJ2 in Eq.Q]. Therefore, the theoretical 
results shown in Figs. [5HT51 depend essentially on the 
model used for describing the production mechanism. 
As discussed above, we used the assumption that the 
7P — > TT^'K^p reaction proceeds exclusively through the 
transition Att — )> ttttN. Agreement with the measure- 
ments would support this asumption. 

In Figs. [8] and [HI we show our experimental results for 
the 7p — )■ 7r°7r°p differential cross sections as a function of 
the production angle of the outgoing tt'^ and proton in the 
center of mass frame. FiguresfTOl and [TTlshow our differ- 
ential cross sections as a function of the invariant mass 
squared rri^ {n'^ tt'^ ) and m?{'K^p). These results are ob- 
tained for the same energies that were used for the angu- 
lar distribution, shown in Figs. [D and [Sj The model pre- 
dictions are shown in these figures up to — 775 MeV. 
They are in reasonable agreement with the experimental 
results, especially concerning the shape of the angular 
distributions. In particular, the calculations reproduce 
not only the convexity and its sign, which changes with 
energy, but also the forward-backward asymmetry. This 
asymmetry is mostly determined by the interference of 
the 3/2~ wave with positive-parity waves (in our case 
1/2+ and 3/2+). As discussed above, the smallness of 
this forward-backward asymmetry indicates the weakness 



of this interference. This is also related to the small val- 
ues of the moments M^sm, shown in Fig.[6l 

The model predictions for the invariant-mass distri- 
butions, shown in Figs. [10] and lll[ are not so impres- 
sive. The poorer level of agreement could be partially 
explained by pion rescattering in the final state, which 
was neglected in our model. According to Refs. |4lll4a|. 
the pion loops in the ■n^ir^ channel can lead to a signif- 
icant enhancement of the cross section at low energies. 
This is primarily because of a large yield of tt+tt^ pairs, 
which in turn can rescatter into neutral pions. As known, 
the interaction between pions in the state = 0+, / = 
is attractive. The corresponding phase shift reaches 7r/2 
close to A'U^ = 900 MeV (see Ref. [is"]). Therefore, it 
is reasonable to expect that the inclusion of this effect 
will shift the m(7r''7r'') spectrum to higher masses. An- 
other possible reason for the poor agreement is that the 
fraction of the wave = 1/2+ in our model is slightly 
overestimated. If the tttt system does not resonate (or the 
TTTT resonance is wide, like /o(600)), then the shape of the 
Dalitz plot (M^^, M^p) is totally determined by the spin- 
parity of a given partial wave (see the corresponding 
discussion for tt^t; photoproduction in Ref. [i3|)- To il- 
lustrate this statement, the contributions of the individ- 
ual states to the tttt spectrum are shown in Fig.[T2] for 
E^f — 425 and 775 MeV. The predicted enhancement of 
the cross section at the boundaries of the kinematical re- 
gion, which is typical for the contribution from the state 
— 1/2+, is not exhibited by the experimental data. 
In this respect, our experimental results prefer a 1/2+ 
fraction that is even smaller than predicted by the fit of 
the moments Wlm- 

In Figs. [5| to [TTJ our experimental results are also 
compared to the predictions of the Bonn-Gatchina model 
[HI, [11]. The approach of Refs. [HI [ll is based on the 
event-by-event likelihood fit that allows one to take accu- 
rately into account the correlations between the different 
reaction channels, for example, ttA and aN . As one can 
see, the model from Refs. [ll|, [11] describes our experi- 
mental results quite well above 550 MeV. Only at low E^ 
does it overestimate the measured cross sections. 

In Fig.[T21 the prediction of our model for the Act = 
(T3/2 ^ (^1/2 helicity asymmetry is compared to the ^p — > 
tt'^'k^P experimental data from Ref. 0]. The experimental 
data, measured as the difference between the total cross 
sections with the initial 7p-system helicity 3/2 and 1/2, 
indicate the dominance of the A = 3/2 component over 
A = 1/2 in the energy region W = 1400 - 1500 MeV, 
excluding any large contribution from the J = 1/2 
waves. Our model reproduces the general trend of the 
data, which shows 0^/2 dominance. The isobar model of 
Ref. , in which the Roper resonance was rather impor- 
tant in the region E^ = 500 — 600 MeV, predicts negative 
values for Act (shown by the dashed line in the same fig- 
ure). In this energy region, this is in contradiction with 
the experimental data. 
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FIG. 13: (Color online) Helicity asymmetry Act — CT3/2 — o"i/2 
for 7p — >■ 7r'^7r''p. The experimental data from Ref are 
shown by open circles. The prediction is shown by the solid 
line for our model, and by the dashed line for the isobar model 
ofRef. il^. 



VI. SUMMARY AND CONCLUSIONS 

The 7p — >■ TT^TT^p reaction has been measured at the 
tagged-photon facility of the Mainz Microtron MAMI-C 
using the Crystal Ball and TAPS spectrometers. The ex- 
perimental results, obtained from the production thresh- 
old up to a photon energy of 1.4 GeV, include the total 
cross sections, various differential cross sections, and spe- 
cific angular distributions. The moments Wlm obtained 
from these angular distributions were used to study the 
importance of different partial waves in double- tt*' pho- 
toproduction at energies below 1)13(1520) (a region that 
has not been so far fully understood theoretically). The 



reliability of our model was checked by the comparison of 
its predictions with the 7p — > tt'^tt'^p experimental results 
for other observables. 

Our analysis of the energy dependence of Wlm showed 
that a large contribution from the J = 3/2 waves is nec- 
essary not only in the region of Z3i3(1520) but also at 
energies below. According to our results, these waves 
seem to be responsible for an almost linear rise of the 
7P — tt'^tt'^p total cross section in the region = 
450 — 725 MeV. Isobar models with the dominant contri- 
bution from 1)13(1520) and a moderate role for the Roper 
resonance cannot explain such features in double- tt'^ pho- 
toproduction. Whether these features are the reflection 
of a large — 3/2+ fraction of tt+tt" — >■ tt^tt" rescat- 
tering, or are a consequence of the strong 1)33(1700) ex- 
citation, found in Ref. [l^ . requires further experimental 
and theoretical studies. 
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